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^! 1 Abstract 
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■ Small neutrino masses can arise in some grand unified models or superstring 

theories. We consider a model with an enhanced fermion sector containing 
Dirac neutral heavy leptons. The dependence on the mass and mixing param- 
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• eters of these new fermions is investigated for several measurable quantities. 

We study the flavor-conserving leptonic decays of the Z boson and universal- 

X 

■ ity breaking in these decays. We also consider the W boson mass dependence 

on neutral heavy lepton parameters. 
PACS numbers(s): 14.60.-z, 12.15.Ff, 13.38.Dg 
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I. INTRODUCTION 



Experimental evidence suggests that neutrino masses are very small, if not zero. How- 
ever, the way the standard model of electroweak interactions (SM) accommodates massless 
neutrinos, by the absence of right-handed neutrino fields, is considered unnatural. Inter- 
esting solutions to this problem have been suggested in the low energy limit of superstring 
theories M, and in grand unified theories (GUT's) ||. Naturally small neutrino masses may 
even be accommodated within the framework of the SU(2)l x U(1)y symmetry. Either the 
fermion content of the SM or the Higgs sector may be extended. In this work, we concentrate 
on the former option. 

We consider an SU(2)l x U{1)y based, superstring-inspired model with an extended 
fermion sector including neutral heavy leptons (NHL's) The NHL's are Dirac particles 

and B-L (baryon minus lepton number) conservation is imposed as an unbroken symmetry. 
This is in contrast to see-saw models [Q,[5] wherein both light neutrinos and NHL's are 
Majorana particles and B-L is broken. 

The model considered allows for the possibilities of lepton-flavor violation, universality 
violation, and CP violation. Our primary interest here is in obtaining potential constraints 
on the model from current experimental data on Z leptonic decay widths. We focus on the 
direct contribution of NHL's to flavor-conserving leptonic Z decays via one-loop diagrams. 
In addition to the Z partial widths for individual lepton flavors, we also calculate a measure 
of the leptonic universality violation, as defined in Ref. 6. Further, we find that we must 
take into account the impact of NHL's on the mass of the W boson since Mw is an input 
parameter within the renormalization scheme adopted here. Hence, we display also the 
dependence of Mw on the parameters of our model. 

This paper is organized as follows. In Sec. U, we describe the model. Existing experi- 
mental constraints on the parameters of the model, namely the masses and mixings of the 



neutral leptons, are reviewed in Sec. |T|. Our one-loop calculation of the Z leptonic decay, 
Z — > l + l~ , is presented in Sec. [TV]. This section also contains a discussion of our renormaliza- 
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tion scheme, including the consideration of M\y- Many of the detailed results are relegated 
to an Appendix. In Sec. [V], we present our results on the Z leptonic widths, the universality 
violating measure, and the W mass. We summarize and draw our conclusions in the final 
section. 



II. DESCRIPTION OF THE MODEL 



Originally, only broken B-L symmetry and Majorana neutrinos were thought of as pro- 
viding an understanding of the smallness of neutrino masses. However, as discussed in Ref. 
4, superstring inspired models can have small neutrino masses (in fact, zero) even if B-L 
symmetry is unbroken and NHL's are Dirac particles. In these models, the SM particle 
content is extended by two new neutrino fields, N R (0, 0) and Sl(0, 0), per family; the zeros 
indicate SU{2) L x U(l)y quantum numbers. Imposing total lepton number conservation 
leads to the mass matrix 



-^- L N L S L ) 



DO 
D T M T 
MO 



N, 



\ Sr J 

where = (z/£, v T L ) and vr = \CPT)v^. D and M are 3x3 mass matrices. The 
diagonalization of the mass matrix yields three massless neutrinos (z/j) along with three 
Dirac NHL's (N a ) of mass ~ M. Note that this implies there are no time dependent 
neutrino oscillations and no neutrinoless double beta decays. The weak eigenstates vl are 
mostly massless neutrinos with a small mixing (~ D/M) of NHL's. The NHL mixing in this 
model is not restricted by small neutrino masses (as is often the case with see-saw models 
where both the mixing and the masses of light neutrinos are sensitive to the D/M ratio), and 
hence rates for all interesting phenomena can be large [f§,[f]|§ . This model is thus attractive 
not only conceptually, but also practically. 

The weak interaction eigenstates Ul are related to six mass eigenstates n a via a 3 x 6 
mixing matrix K with components Ki na ; I = e, /x, r and n a = u±, V2, V3, iV 4 , N 5 ,N e . 



(1) 



K eVl K eV2 K eV3 K eN4 K eNr> K eN& 
Ki na = K^ V2 K^ Na K^. = (K L , K H ). (2) 

^ K TU1 K TU , 2 K TL , 3 K t n 4 K t n 6 K t n 6 

After rotating away redundant degrees of freedom from K, we are left with 3 2 angles and 
(3 — l) 2 phases. This allows for possible lepton-flavor violation, universality violation and 
CP violation. 

The mixing factor which typically governs flavor- conserving processes, is given by: 



n a =N 4 ,N 5 ,N 6 

and the flavor-violating mixing factor IJbmix is defined as: 

lahmix = ■^■l a n a ^-l b n a ^ ^ a f> = e > 7^ ^b- 

n a =N 4 ,N 5 ,N 6 



(3) 



(4) 



Further, an important inequality holds: 



1/ / 2 I < / / / / 



a^b. 



(5) 



This implies that one might observe nonstandard effects in flavor-conserving processes even 
if they are absent in flavor-violating processes. 

For reference, the charged current Lagrangian is given by 

£ C c = -^j=gW IJ ' hiii 1 - l5)Ki na n a ] n a = u u v 2 , ^3, N 4 , N 5 , N 6 (6) 

and the neutral current Lagrangian as (the ZNN part is obtained by analogy) 



C nc = -2-Z" £ ^(kIKh)^! - 75 )iV a , 



4c 



W i=l,2,3;a=4,5,6 

where cw = cos9w, ®w being the Weinberg angle. 



(7) 



III. REVIEW OF EXISTING CONSTRAINTS ON NEUTRAL HEAVY LEPTONS 



Constraints on neutral heavy lepton masses and mixings come from three different 
sources. First, there is the possibility of direct production of NHL's. For instance, if an 
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NHL is light enough, it could be produced in some decays, e.g. Z — > iV a + u, and subse- 
quently decay itself. The rate for Z decays into an NHL and a light neutrino has been given 
previously as 

M 2 M 2 

T(Z ^N a + u) = a mix (l ~ -j^r)(l + + (8) 

where 

dmix = X! 1-^iiVal 2 - (9) 

The subsequent NHL decay rate (for M N < M w ) is then given by 

IV = a mix {^f^Y^ (10) 

where T M is the muon decay rate and is the effective number of decay channels available 
to the NHL ||. Given the absence of experimental evidence for such direct production, we 
will consider only NHLs with mass greater than the Z mass. 

Secondly, there are constraints on NHL mixing parameters from a variety of low energy 
experiments and from experiments at the CERN Large Electron Positron Collider I (LEP I). 
Due to unitarity properties of the mixing matrix K, a nonzero NHL mixing slightly reduces 
the couplings of light neutrinos from their standard model values, thus affecting rates for 
nuclear (3 decays, r and ir decays, and for Z decays. The following upper limits are consistent 
with experiment [[K]] 

ee mix < 0.0071 

Wmix < 0.0014 (11) 
TT mix < 0.033 or < 0.024 including LEP I 

These limits are model independent and also independent of the NHL mass. They arise from 
a global analysis of results including lepton universality experiments, Cabibbo-Kobayashi- 
Maskawa (CKM) unitarity tests, W mass measurements and results from LEP I experiments. 
Note that the LEP constraints presented above do not include NHL loop effects but, rather, 



only coupling constant modifications due to mixing. We consider NHL loop effects in this 
work. 

Finally, the NHL masses and mixings can be constrained via their contribution in loops 
to various processes. Such constraints are NHL mass dependent. Flavor-violating processes 
have previously been studied at low energies; these include \i — > ej, \i — > 3e Q,|,|TT|] and 



flavor violating decays of the r H,^- F° r instance, in the context of the model considered 



here, with mass degenerate NHLs, the /i — > branching ratio is @|| 

BR(fi -> ej) = ■^r-\efimi X \ 2 \F^(x)\ 2 (12) 

OZ7T 

M 2 

where x = j^f, and F 7 {x) is an NHL mass dependent form factor. For NHL masses Mm > 
500 GeV, which we will ultimately consider, F 7 (x) — > —2. Given the current experimental 
limit on the /z — > branching ratio (< 4.9 x 10~ n ) [[L|], this yields an upper limit on the 
mixing of 

|e/w| < 0.00024. (13) 

By combining the constraints obtained from the global analysis (Eq. (11)) with the inequality 
relations of Eq. (^) one obtains the following upper limits on the mixing factors 









| (^fJ'mix 


< 0.0032 








| f^Tmix 


< 0.0068 








| &T~mix 


< 0.015 


For the mixing 


JS [lT m i x 


and 


GTmix, these are the stron 


flavor-violating 


I leptonic Z 


decays, Z — > 


e/i, er, fir, 



(14) 



work, we consider the flavor-conserving decays Z — * ee, /i/i, rr. 
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IV. CALCULATION OF THE ONE-LOOP LEVEL CONTRIBUTION OF 
NEUTRAL HEAVY LEPTONS TO LEPTONIC FLAVOR-CONSERVING Z 

DECAYS 



As noted previously, the limits on mixing parameters extracted from LEP I observables 
1(J do not include NHL's in the one- loop diagrams; only mixing factor modifications are 



made to the tree level results. We consider here the direct contribution to flavour conserv- 
ing leptonic Z decay of NHL's via one- loop diagrams. The importance of studying these 
processes is enhanced by the implication of Eq. (H); one may observe flavor- conserving (but 
universality breaking) Z decays even in the absence of flavour violating processes. 

For Z leptonic decay, NHLs contribute directly to Z oblique corrections, as shown in 
Fig. |l|, to lepton wave function renormalizations and to vertex corrections, as in Fig. |2|. 
These one- loop contributions of NHL's can be incorporated into the framework of the full 
standard model one-loop electroweak corrections. The one-loop corrected leptonic width 
can be parametrized as |15] 

rz= ^f° (i + s QED ), (15) 
i + n z (M|) 

where the tree level leptonic width of the Z boson is given by 

r Q = ^M z (v} + a}), (16) 

Vf and a/ being respectively the vector and axial vector couplings of charged leptons to Z. 
The one-loop electroweak corrections include 5Tz, which represent vertex loops, and flz and 
Sqed which represent the Z-oblique corrections and QED corrections respectively. 

Our calculation is done within the framework of an on-shell renormalization scheme as 
detailed in Ref. 15. All the SM one-loop diagrams were calculated using standard routines 
from the CERN electroweak library |15] modified by appropriate mixing factors. The vertex 
parameter 5T z actually includes also fermion wave function renormalization and counterterm 
contributions in the scheme we adopt. This on-shell renormalization scheme takes a, Mz 
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and Mw as input parameters. However, the direct measurement of Mw is not yet precise 
enough for its use as an input parameter; hence, it is replaced by via the one-loop relation 

^VK s iy = (1 — Ar) ~~ 2 e6mi:c — ij^^raix), (17) 

where Ar is an equivalent of the SM quantity Ar SM , and sw = sin^jy. As a result, we 
also have to consider muon decay loops with NHL's. The types of corrections involved in 
muon decay are pictured schematically in Fig. 3. They include boxes and vertices with 
NHL's as well as the lepton wave function renormalizations and W oblique corrections. The 
NHL diagrams contributing to the W oblique corrections are shown in Fig. 1 along with 
the Z oblique corrections. These corrections all feed into the Z leptonic decay calculation 
indirectly via the dependence of My/ on the parameter Ar and the overall factor modified 
by mixings. 

Referring to Eq. (15) now, the QED corrections, Sqed, are not modified from the SM. 
The factor (1 + n^(M|)) represents the wave-function renormalization of the Z boson. It 
depends on all the unrenormalized propagator corrections (E^ , T*w , 'H-yZ, E 7 ) fTSfl . Of these, 
X 7 ^ and S 7 are not modified within our model while Hz, Hw both contain nonstandard 
terms. Those nonstandard terms, denoted as (s) and H w N (s), are given by Eqs. (A2) 
and (A4) respectively, in the Appendix. They consist of Mn dependent terms representing 
the direct contribution from NHL's in loops in Fig. 1 and of SM terms modified by mixing 
factors (the indirect effect of NHL's reducing the mixings of light neutrinos through the 
unitary matrix K). 

The vertex parameter 5T z includes 7 — Z mixing, external fermion wave function renor- 
malization, and counterterm contributions in addition to the vertex loops involving NHL's. 
Those fermion self energy and Zff vertex loops which contain NHLs are shown in Figs. 
2a-j. The individual contributions of Figs. 2a-j are given in Eq. (A6) in the Appendix. 

It is characteristic that NHL's in loops generally do not decouple (violation of the 
Appelquist-Carazzone theorem) ||16|1 ; rather they often show a quadratic mass dependence 
~ -gf-- This is a common feature for theories based on the spontaneous symmetry breaking 



mechanism. We are already familiar with a similar result for the top quark mass in SM loops 
fTTfl . Indirect bounds on the top quark mass arise from the Z and W polarization diagrams 
and Zbb vertex loop since these corrections come in with the amplitude 



2 2 

M <~^M^W (18) 



While the top quark mixes with the full strength (\V tb \ ~ 1), the NHL mixings are limited 
by Eqs. (11) and (14). As a result, we find sensitivity only for ~ 10 m t . Thus we will 
only present numerical results for NHL masses greater than about 500 GeV. 

Given the nondecoupling feature, the vertex correction 5T z is dominated for large NHL 
masses by the diagrams (in decreasing order of importance) 2j, 2f and 2e, while 2a-b, 2c-d, 
2g, 2h and 2i are negligible (the largest unrenormalized contribution comes from 2g and 
2c-d; however, diagrams 2c-d enter the renormalized vertex correction STz as a part of a 
counterterm that cancels out the large amplitude of the graph 2g). 

To illustrate how flz and Ar depend on M N (for M N large), we separate M N dependent 
terms as Hm n , ^m n and find in the limit X~ = j^- — > (large NHL mass) 



Ar Mw : 



TVS 



W 



16 24 , mf 

i=e,/i,T ' 



a 

71 



where 



r 2 - <? 2 1 M 2 

Ihh W ? W X + (Ill-1)— -t y In^f 

i0S W Z4 "V l=e,ti,r m l 



|2 



(19) 



l=e,/j,,T 

Ill = 1 - 2 (ee mix + n^ mix + TT mix ) + l HH (20) 

With the NHL mass of the order of several TeV, one has to worry about the perturbative 
unitarity bound. A good way to demonstrate this is to bring about the Higgs analogy. The 
width of the Higgs boson of mass m# is given by 

T H = (21) 
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which can be compared with the width of the NHL. For NHL mass Mjy 3> M\v. M/. Mn- 
the width is fl~8|j : 



a 



N 



M%a mix . (22) 



Demanding that Th < \M H we get the well known bound on the Higgs mass M H < 1 TeV 
Similarly, demanding < \Mn for the NHL, with the current constraints on mixings, (Eq. 
(11)), one obtains M N < 3.5 TeV. 

V. RESULTS 

In this Section, we present our numerical results. As input parameters, we used Mz = 
91.173 GeV, M H = 200 GeV, a' 1 = 137.036 and A = ^g- = 37.281 GeV. We have 
assumed degenerate masses for the three NHL's and present results for the NHL mass range 
0.5 TeV < Mat < 5 TeV, as motivated by the non-decoupling and perturbative unitarity 
arguments given in the last Section. We have also imposed restrictions on the mixing 
parameters. We assume that ee m i X and fi/i m i x are very small relative to TT mix . The model 
and NHL mass independent limits quoted in Eq. (11) are more stringent for e and \i than 
for r. In addition, our assumption is also partially supported by the smallness of efi mix , as 
determined from /i — > cy, in combination with the inequality Eq. (5). This neglect of ee m i X 
and mimix proves useful practically in that many of the muon decay loops (boxes and vertex 
corrections, but not W oblique correction) are eliminated as a result. 

The Z leptonic width is given as a function of NHL mass in Figs. 4a, b. In Fig. 4a, 
we have fixed the mixing TT mix = 0.033. The width for Z decay to e + e~ is shown for a top 
quark mass of 174 GeV. The Z decay rate into r + r _ is shown for three values of the top 
quark mass, 150, 174 and 200 GeV. The dashed lines represent the la variation about the 



current experimental result for the average Z-leptonic width of T/ = 83.96 ± 0.18 MeV [|T9 
In Fig. 4b, we fix the top quark mass at 174 GeV and show the Z width to t + t~ for three 
values of the mixing parameter, TT mix = 0.02, 0.033, 0.07. We also present results for the 
universality breaking ratio defined as || 
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u. 



br 



T(Z -> T+T-) - T(Z -> e + e~) 



(23) 



r(Z -> r+r-) + r(Z -> e+e" 

This is shown in Fig. 5 as a function of M/y, again with m 4 = 174 GeV, and the mixing 
parameter varied about TT mix = 0.033. The la experimental limit on U& is indicated as 
the dashed line. Note that the most recently reported values of Z widths into individual 
lepton flavors (|TUJ have T(Z — > r + r _ ) > T(Z — > e + e~), Y{Z — > as opposed to the 

last round of results |L3| . 

Finally, we present the NHL mass dependence of the W mass in Figs. 6a, b. The top 
quark mass is varied in Fig. 6a, while the mixing is held constant at TT mix = 0.033. In Fig. 
6b, the mixing is varied about rr m j X = 0.033 for a fixed top quark mass m t = 174 GeV. 



VI. DISCUSSION AND CONCLUSIONS 

Our primary consideration here has been the inclusion of neutral heavy leptons in the 
calculation of the flavor-conserving Z decays to charged leptons at one-loop level. The 
dependence of the Z leptonic widths on the NHL mass, Mat, and on the mixing parameter 
TT mix which we retain, was given in Figs. 4 a,b. We see for the experimentally allowed upper 
limit of TT m i X = 0.033, and assuming a top quark mass mt = 174 GeV, the Z decay width 
to t leptons is sensitive at the present 2a level to NHL masses larger than about 2.5 TeV. 
The top mass dependence is also shown in that Figure. The sensitivity to and mt arises 
since these heavy fermions generally do not decouple from the one-loop diagrams. Fig. 4 
b indicates how the Z width dependence on Mjy varies with the mixing parameter. Apart 
from this comparison of each leptonic width prediction with experiment we can also exploit 
the flavor universality violation which takes place in the model. The universality breaking 
ratio, Ubr, defined in Sec. 4, is sensitive to NHL masses above approximately 3.5 TeV at the 
la level, assuming TT mix = 0.033. 

The W boson mass also exhibits some sensitivity to NHL parameters arising from the 
mixing factor modifications and the presence of one-loop diagrams containing neutral heavy 
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leptons, as described in Sec. 3. From Figs. 5a, b we see that the W mass, currently measured 
as Myy = 80.23 ± 0.18 GeV |[20]| , is sensitive at the la level to NHL masses greater than 
about 3.5 TeV, again assuming TT m i X = 0.033 and mt = 174 GeV. The experimental error 
on Mw might be expected to come down to about 0.05 GeV once LEP II measures W pair 
production |?T . 



We have considered a model containing isosinglet neutral heavy leptons which can ac- 
commodate various phenomena beyond the standard model, such as lepton flavor-violation, 
CP violation and lepton universality violation. We have presented the dependence on the 
mass and mixing parameters of this model for Z decays to charged leptons and for the W 
boson mass. Because the NHL mass and mixing dependence is different for the Z decay 
width and the W boson mass, they provide somewhat complementary information on these 
parameters. 

Current data from LEP I on Z leptonic widths and the present Collider Detector at 
Fermilab (CDF) and DO Collaboration measurements of My/ are sensitive to NHL masses 
greater than about 2.5 - 3.5 TeV. With the accumulation of about 60 (pb)^ 1 at LEP I in 
1994 and the prospect of the very precise W mass measurement at LEP II, these sensitiv- 
ities will certainly be improved considerably. Thus the Z partial width to leptons and W 
mass measurements can provide, along with the other observables discussed in Sec. 3, a 
consistency check on the possible existence of isosinglet neutral heavy leptons. 

ACKNOWLEDGMENTS 

This work was funded in part by the Natural Sciences and Engineering Research Council 
of Canada. We would like to acknowledge useful communications with M. C. Gonzalez- 
Garcia and A. Pilaftsis. 



12 



APPENDIX: 



In this Appendix, we present the parts of our calculation which are necessary to obtain 
the NHL dependent modifications of the one-loop Z decay width to charged leptons. We 
refer the reader to Ref. 15 for further details. We define the following symbols and mixing 
factors: 



X 



2 3 
A = 7 + In An + - In /i 2 , 

6 2 

A m = - - 7 + ln47T-ln^, 

e fx 2 

I2 \l&mix\ \lfJ<mix\ l^mix] i 
^3 ttmix I2 1 
I4 — 1 — lllmix + h, 
IcH &&mix l^l^mix TTmixi 

Ihh = E ^2, 

l=e,fi,r 

Ilh = %{Ich — Ihh), 

l LL = 1 - 21ch + Ihh- (Al) 



First we deal with the direct contributions of NHL's to the Z boson wave function renormal- 
ization, as parametrized by n^(Mf ). fLz(Mf) depends on two quantities which are modified 
from the SM by the inclusion of NHL's. The neutral lepton part of the unrenormalized self 
energy of the Z boson, which consists of the SM massless v loop, the mixed uN loop and 
the NN loop, is given by 

Z u /(s) = -^\ S -Ill E [A"* + 2 - ln(~ - ic) - £ 
8tts 2 w c 2 w {3 l= ^ r [ m 2 3 J 

+ Ilh { S --Mk ) + l s _Mk + F{s , 0, Mn) ( S - - M _ <)1 
L v 3 2 ; 9 6 v A 3 6 6s '. 
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+ Ihh 

where the function F(s, mi, m-i) is given by 



A M -(^ - Ml) - S - + ±F(s, M N , M N )(s - M 



. . to? + mo , m i 

t is, mi, m 2 ) — — H 5 5- m — 

to| — to 2 to 2 



f 1 — x(s + TO? — TOo) + to? — it 

— / oxln . 

JO TOxTO 2 



(A2) 



(A3) 



The leptonic part of the unrenormalized self energy of the W boson is given by: 



a 



+(*- 

+ E h 

l=e,n,T 



I 


E (i 






TO; 2 


4 

TO; 


2 


2s 








[s - — to 2 ) A mi 



)F( a ,0,m,) + | a -^ 
5, „ to 2 A m < 



[S — -TO; 



+ ( S -^-(M1!)^M,,to,) 



+(*- 



M 2 N + mf 



)(1 



2s 

M^+to? M 



'v 



In 



iV V 



TO; 



TO; 3. 



(A4) 



Next, we consider the contributions of NHL's to the parameter 5T Z . The direct contri- 
bution of NHL's in the triangle diagrams of Figs. 2e-j is given by the sum of amplitudes: 

a 



CM ( 

M = +iee^(l - 75)— \hMzsw + hM Z NN* + hM z ** + (1 - k)M 

47T L 



ZWWv 



-Ii-M-zwwn + hM-zNvw + h-MzvNw + h-M-zww + h-M-zNNw}, (A5) 



with (jS- terms neglected) : 



•M.Z&W — + 



2s\yCw 



XC (M W ,M N ,M W ), 



•nA-ZNN® — -T 



8s^c w 
2s^ 2c^ 



^ 2 Co(M^,M w ,M^), 



X[Cit(M w ,M N ,M w ) + ^A], 
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M ZW wa = ^M* z [-C n {M Wi M ai M w ) 

-C 23 (M W , M a , M w ) - C (M Wl M a , M w )\ 
+4CH n (M w , M a , M w ) - ~ + A} , 
Mzabw = \2M 2 z [(C 23 (M a , M w , M b ) 

oS w Cw L 

+C n (M a , M w , M h )\ + 2 - 4C|f(M , M 6 ) - A} 

where a,b run over N, v ; M v = 0. (A6) 



Here Aiz$w is the sum of equal contributions from diagrams 2h and 2i. Diagram 2f comes in 
both with massless z/'s and NHL's. Diagram 2e comes in with four combinations of neutral 
lepton types. Thus Aizwwu and M.zww are standard model results but come into the sum 
(Eq. (A5)) with NHL mixing factor coefficients. Our results in Eq. (A5) are written in 



terms of the 't Hooft-Veltman integrals \T2\. Our conventions are given below, with finite 
parts indicated by the superscript. 
The function C is defined as: 

Co(mi,m 2 , 777,3) = C (pi,p 2 ; mi, m 2 , m 3 ) = Co m (m 1 ,m 2 ,m 3 ) 

d n q 1 



m 2 D' 

where D = (q 2 — m\ + ie)[(q — Pi) 2 — m\ + ie] 

x[(q-Pi -P2) 2 -ml + ie]. (A7) 



The functions C24, C23, C\\ are defined by: 

/d n q q^ 
~2~p = ~V\liP\\ — p 2(J C\ 2 , 



a 



d n q q^q v 



^ J in 2 D 

= V\ i iV\vC 2 \ + V2 ll ,V2vC 22 + {j>\iiV2v + V\vV2^C 23 - g^ u C 2i . (A8) 
The functions C24, C 23 , C\\ are reduced (in the limit p\ = p\ = mf <C (pi + p 2 ) 2 = M|) to: 
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1 2 t- 
C 2 4(m 1 ,m 2 ,m 3 ) = ^(j— - - 7 - bi7r) + C%r(m 1 ,m 2 ,m 3 ), 

C£l n (m 1 ,m 2 ,m 3 ) = [mfC (m 1 ,m 2 ,m 3 ) + / 1 C u (m 1 ,m 2 ,m 3 ) + 

fin 1 1 

5f n (pi +p 2 ;mi,m 3 )](--) + - 



with fi=m 1 —m 



2' 



fin 1 

C 11 (m 1 ,m 2 ,m 3 ) = C(l n (m 1 ,m 2 ,m 3 ) = -j^ 2 -[f 2 C (m 1 ,m 2 ,m 3 ) - 

z 

Bo m (Pi + P2; m 1 ,m 3 ) + B^ n (pi, m 1 , m 2 )\ 



with / 2 = M| + Wlj — 77^3, 
C 23 (mi, m 2 , m 3 ) = C^ n (mi, m 2 , m 3 ) = --^-[.B^ m (pi + p 2 ; mi, m 3 ) + B f Q m (p 2 ] m 2 , m 3 ) 



2 

+/iCn(mi,m 2 ,m 3 )] + ^"(mi, m 2 , m 3 )— j. (A9) 



The functions i?o, -£>i are defined as: 



B (p;m 1: m 2 ) 



in 2 (q 2 — m 2 + ie) [(q — p) 2 — m 2 . + it] 
2 

= (-: 7-lnvr) + J B^ n (p;m 1 ,m 2 ), 

4 — n 



■>Q m (p;mi,m 2 ) = — dxln\p 2 x 2 + m\ — (p 2 + m\ — m%)x], 
■In 

1 2 

5i(p;mi,m 2 ) = --(- 7 - lnvr) + B( m (p;m 1 ,m 2 ), 

2 4 — n 

?{ m (p; mi, m 2 ) = / drr ln[p 2 x 2 + m 2 — (p 2 + m 2 — m\)x]x. (A10) 

JO 



The remaining diagrams of Fig. 2 contribute NHL dependent terms to the charged lepton 
wave function renormalization. The contributions to the lepton self energy of Figs. 2a, b 
and Figs. 2c, d are given as Y> w and £*, respectively. 

- iZ w = -^Hi[2 - 2 ln4vr - 2 ln/i 2 + 2 7 - - + Af(X)}p al a (l - 7s ), 
Slixsiw e 

ia 2 

= -^r l i X [- + ln4vr + ln/i 2 - 7 - 2/(*)K7 a (l - 75 ), 

where f(X) = , — NO In X + — * + 1 - + -InM 2 . (All) 

Jy ! 2{X-l) 2 1 — X 4(1 — X) 2 w V ; 
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All these contributions, Eq. (A5) and (All), along with their corresponding counterterms 
modify the Zl + l~ vertex. In addition, the vertex is modified via 7 — Z mixing. However, 
the relevant term IT 7Z (M§ ) depends on NHL's only through the Z and W self energies so 
we need no other results to determine 5Tz- 

The results presented in this Appendix are also sufficient to derive the W oblique cor- 
rections of Ar, which we have considered. 
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FIGURES 



FIG. 1. Diagrams for oblique corrections due to neutral heavy leptons N. 

FIG. 2. Diagrams for one-loop vertex correction to flavor-conserving leptonic Z decays due to 
neutral heavy leptons N. 

FIG. 3. Schematic representation of the one-loop muon decay diagrams with neutral heavy 
leptons 

FIG. 4. Z leptonic width as a function of Mn for (a) fixed mixing parameter and different 
values of m< (b) fixed m t and different values of the mixing parameter. The dashed lines represent 
Id band about the current experimental value Ti = 83.96 ± 0.18 MeV. 

FIG. 5. Universality breaking parameter Ub r as a function of for fixed mt and different 
values of the mixing parameter. The dashed line represents la experimental limit (< 0.005). 

FIG. 6. W mass as a function of Mn for (a) fixed mixing parameter and different values of mt 
(b) fixed mt and different values of the mixing parameter. The dashed lines represents la band 
about the current experimental value M\v = 80.23 ± 0.18 GeV. 
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